1 



STELLAR STRUCTURE AND 
EVOLUTION: DEDUCTIONS FROM 

o HIPPARCOS 

o ■ 

■ 

(N : 

5h ■ Yveline Lebreton 

1 c3 ■ 

Observatoire de Paris, DASGAL-UMR CNRS 8633, Place J. Janssen, 92195 Meudon, France 
£^ , e-mail: Yveline.Lebreton@obspm.fr 



in 

' KEYWORDS: Stars: physical processes, distances, H-R diagram, ages, abundances, Hipparcos 

o : 
o . 

O • ABSTRACT: 

O,; 

q ■ evolution of stars has become possible as a result of enormous progress made in 
H ' 

| the acquisition of high-quality observational and experimental data and of new 



During the last decade, the understanding of fine features of the structure and 



developments and refinements in the theoretical description of stellar plasmas. 
The confrontation of high-quality observations with sophisticated stellar models 



^ | has allowed many aspects of the theory to be validated, and several character- 
istics of stars relevant to Galactic evolution and cosmology to be inferred. This 
paper is a review of the results of recent studies undertaken in the context of 
the Hipparcos mission, taking benefit of the high-quality astrometric data it has 
provided. Successes are discussed, as well as the problems that have arisen and 
suggestions proposed to solve them. Future observational and theoretical devel- 
opments expected and required in the field are also presented. 
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1 INTRODUCTION 

Stars are the main constituents of the observable Universe. The temperatures 
and pressures deep in their interiors are out of reach for the observer, while the 
description of stellar plasmas requires extensive knowledge in various domains 
of modern physics such as nuclear and particle physics, atomic and molecular 
physics, thermo- and hydrodynamics, physics of the radiation and of its interac- 
tion with matter, and radiative transfer. The development of numerical codes to 
calculate models of stellar structure and evolution began more than forty years 
ago with the pioneering works of Schwarzschild (1958) and Henyey et al (1959). 
These programs have allowed at least the qualitative study and understanding 
of numerous physical processes that intervene during the various stages of stellar 
formation and evolution. 
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During the last two decades, observational data of increasingly high accuracy 
have been obtained as a result of 1) the coming of modern ground-based or space 
telescopes equipped with high-quality instrumentation and with detectors giving 
access to almost any possible range of wavelengths and 2) the elaboration of var- 
ious sophisticated techniques of data reduction. Ground-based astrometry has 
progressed, while space astrometry was initiated with Hipparcos. In the mean- 
time, CCD detectors on large telescopes opened the era of high-resolution, high 
signal-to-noise ratio spectroscopy while multi-color filters were designed for pho- 
tometry. New fields have appeared or are under development, such as helio- and 
asteroseismology or interferometry. On the other hand, stellar models have been 
enriched by a continuously improved physical description of the stellar plasma, 
while the use of increasingly powerful computers has led to a gain in numerical 
accuracy. 

The confrontation of models with observations allows testing and even valida- 
tion of the input physics of the models if numerous observations of high quality 
are available. Fundamental returns are expected in many domains that make use 
of quantitative results of the stellar evolution theory such as stellar, Galactic, and 
extragalactic astrophysics as well as cosmology. Because of their positions, move- 
ments, or interactions with the interstellar medium stars are actors and tracers 
of the dynamical and chemical evolution of the Galaxy. Astrophysicists aim to 
determine their ages and chemical compositions precisely. For example, the firm 
determination of the ages of the oldest stars, halo stars or members of globular 
clusters, is a long-standing objective because it is one of the strongest constraints 
for cosmology. 

Although great progress has been made, a number of observations cannot be 
reproduced by stellar models, which raises many questions regarding both the 
observations and the models. In the last few years, two scientific meetings have 
been explicitly devoted to unsolved problems in stellar structure and evolution 
(Noels et al 1995, Livio 2000). A major point of concern is that of transport 
processes at work in stellar interiors (transport of the chemical elements, angu- 
lar momentum or magnetic fields by microscopic diffusion and/or macroscopic 
motions). Observations show that transport processes are indeed playing a role 
in stellar evolution but many aspects remain unclear (sometimes even unknown) 
and need to be better characterized. Another crucial point concerns the atmo- 
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spheres, which link the stellar interior model to the interstellar medium and are 
the intermediate agent between the star and the observer. Uncertainties and 
inconsistencies in atmospheric descriptions generate errors in the analysis of ob- 
servational data and in model predictions. 

This paper is the third of the series in ARAA dedicated to the results of the 
Hipparcos mission; Kovalevsky (1998) presented the products of the mission and 
the very first astrophysical results obtained immediately after the release of the 
data, while Reid (1999) reviewed the implications of the Hipparcos parallaxes for 
the location of the main sequence (MS) in the Hertzsprung-Russell (H-R) dia- 
gram, the luminosity calibration of primary distance indicators, and the Galactic 
distance scale. Also, van Leeuwen (1997) presented the results of the mission, 
and Baglin (1999) and Lebreton (2000) discussed the impact of Hipparcos data 
on stellar structure and evolution. 

Hipparcos has provided opportunities to study rather large and homogeneous 
samples of stars sharing similar properties, for instance, in terms of their space 
location or chemical composition. I review studies based on Hipparcos observa- 
tions which (1) confirmed several elements of stellar internal structure theory, 
(2) revealed some problems related to the development of stellar models, and (5) 
yielded more precise characteristics of individual stars and clusters. In Sections 
2 and 3, I discuss the recent observational (including Hipparcos) and theoreti- 
cal developments from which new studies could be undertaken. In Section 4, I 
concentrate on the nearest stars, observed with highest precision (A-K disk and 
halo single or binary field stars, and members of open clusters). In Section 5, I 
review recent results on variable stars, globular clusters and white dwarfs based 
on Hipparcos data. The stars considered are mostly of low or intermediate mass, 
and except for white dwarfs, the evolutionary stages cover the main sequence and 
subgiant branch. Throughout this paper, I emphasize that the smaller error bars 
on distances that result from Hipparcos make the uncertainties on the other fun- 
damental stellar parameters more evident; fluxes, effective temperatures, abun- 
dances, gravities, masses and radii have to be improved correspondingly, implying 
in many cases the need for progress in atmospheric description. 
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2 NEW HIGH-ACCURACY OBSERVATIONAL MATERIAL 



This section presents a brief review of Hipparcos results and complementary 
ground-based or space observations which, if combined, provide very homoge- 
neous and precise sets of data. 

2.1 Space Astrometry with Hipparcos 

The Hipparcos satellite designed by the European Space Agency was launched in 
1989. The mission ended in 1993 and was followed by 3 years of data reduction. 
The contents of the Hipparcos Catalogue (Eur. Space Agency 1997) were de- 
scribed by Perryman et al (1997a). The data were released to the astrophysical 
community in June 1997. General information on the mission is given in van 
Leeuwen's (1997) and Kovalevsky's (1998) review papers. 

Stars of various masses, chemical compositions and evolutionary stages located 
either in the Galactic disk or in the halo were observed; this was done systemat- 
ically to a V-magnitude that depends on the galactic latitude and spectral type 
of the star, and more generally, with a limit of V~12.4 mag. The Hipparcos Cat- 
alogue lists positions, proper motions, and trigonometric parallaxes of 117955 
stars as well as the intermediate astrometric data, from which the astrometric 
solutions were derived; this allows alternative solutions for the astrometric pa- 
rameters to be reconstructed according to different hypotheses (see van Leeuwen 
k Evans 1998). 

A total of 12 195 double or multiple systems are resolved (among which about 
25 % were previously classified as single stars), and 8 542 additional stars are sus- 
pected to be non-single. Detailed information on multiple systems, as described 
by Lindegren et al (1997), can be found in The Double and Multiple System 
Annex of the Catalogue. 

The median accuracy on positions and parallaxes (it) is typically ~1 milliarcsec- 
ond (1 mas), whereas precisions on proper motions are about 1 mas per yr. Pre- 
cisions become much higher for bright stars and worsen toward the ecliptic plane 
and for fainter stars. The astrometric accuracy and formal precision of Hipparcos 
data have been investigated by Arenou et al (1995) and Lindegren (1995), and 
discussed by van Leeuwen (1999a): for the Catalogue as a whole, the zero-point 
error on parallaxes is below 0.1 mas and the formal errors are not underestimated 
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by more than 10%. After Hipparcos about 5 200 single stars and 450 double (or 
multiple) stars have parallaxes known with an accuracy cr n /ir better than 5%, 
20 853 stars have cr n /ir lower than 10% and 49 333 stars have a n /n lower than 
20% (Mignard 1997). Martin et al (1997, 1998) and Martin & Mignard (1998) 
determined the masses of 74 astrometric binaries with accuracies in the range 
5-35%. Sderhjelm (1999) obtained masses and improved orbital elements for 205 
visual binaries from a combination of Hipparcos astrometry and ground-based 
observations; among these, 12 (20) systems have mass-errors below 5 (7.5)%. 

The Hipparcos Catalogue also includes detailed and homogeneous photometric 
data for each star, obtained from an average number of 110 observations per star. 
The broad-band Hipparcos (Hp) magnitude corresponding to the specific pass- 
band of the instrument spanning the wavelength interval ~ 350-800 nanometers 
(see Figure 1 in van Leeuwen et al 1997), is provided with a median precision 
of 0.0015 mag for Hp<9 mag. The Johnson V magnitude derived from com- 
bined satellite and ground-based observations is given with a typical accuracy of 
0.01 mag. The star mapper Tycho had passbands close to the Johnson B and V 
bands and provided two-color Bx and Vt magnitudes (accuracies are typically 
0.014 mag and 0.012 mag for stars with Vt < 9). 

Hipparcos provided a detailed variability classification of stars (van Leeuwen et 
al 1997), resulting in 11 597 variable or possibly variable stars. Among these 2 712 
stars are periodic variables (970 new) including 273 Cepheids, 186 RR Lyrae, 108 
S Scuti or SX Phoenicis stars, and 917 eclipsing binaries. 

Hipparcos was planned more than fifteen years ago, and while its development 
proceeded, significant progress was made in the derivation of ground-based paral- 
laxes using CCD detectors. Parallaxes with errors less than 1.4 mas have already 
been obtained for a few tens of stars, and errors are expected to drop to ±0.5 
mas in the years to come (Harris et al 1997, Gatewood et al 1998). In addition, 
the Hubble Space Telescope (HST ) Fine Guidance Sensor observations can pro- 
vide parallaxes down to V~15.8 with errors at the 1 mas level (Benedict et al 
1994, Harrison et al 1999). However, the distances of a rather small number of 
stars will be measured by HST because of the limited observing time available for 
astrometry. The enormous advantage of Hipparcos resides in the large number 
of stars it dealt with, providing homogeneous trigonometric parallaxes that are 
essentially absolute. 
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2.2 Ground-Based Photometry and Spectroscopy 

The fundamental stellar parameters (bolometric magnitude Mboi, effective tem- 
perature T e ff, surface gravity g, and chemical composition) can be determined 
from photometry and/or from detailed spectroscopic analysis. However, the de- 
termination largely relies on model atmospheres and sometimes uses results of 
interior models. Direct masses and radii can be obtained for stars belonging to 
binary or multiple systems. Interferometry combined with distances yields stellar 
diameters giving direct access to T e g, but still for a very limited number of rather 
bright stars which then serve to calibrate other methods. The different methods 
(and related uncertainties) used to determine the fundamental stellar parame- 
ters mainly for A to K Galactic dwarfs and subgiants are briefly discussed, and 
improvements brought by Hipparcos are underlined. 

BOLOMETRIC MAGNITUDES. Integration of UBVRIJHKL photometry gives 
access to the bolometric flux on Earth Fboi, at least for F-G-K stars where most 
energy is emitted in those bands and which are close enough not to be affected 
by interstellar absorption; the (small) residual flux, emitted outside the bands, 
is estimated from model atmospheres. Recently, Alonso et al (1995) applied 
the method to ~100 F-K dwarfs and subdwarfs and obtained bolometric fluxes 
accurate to about 2% and, as a by-product, empirical bolometric corrections for 
MS stars. 

If Fboi and distances are known, Mboi can be derived with no need for bolo- 

i 

metric correction. The accuracy is then 0"M bo i = loge [(2.5 ) 2 + (5^) 2 ] 2 , 
meaning that if -^f ~2% then £JM bol is dominated by the parallax error as soon 
as £^ >1%. In other cases, when the distance is known, Mboi is obtained from 
any apparent magnitude m and its corresponding bolometric correction BC(m), 
derived from empirical calibrations or from model atmospheres. Up to now Hip- 
parcos magnitudes Hp have not been used extensively, despite their excellent 
accuracy (0.0015 mag), because of remaining difficulties in calculating BC(Hp) 
(Cayrel et al 1997a). 

EFFECTIVE TEMPERATURES. The InfraRed Flux Method (IRFM; Blackwell 
et al 1990), applicable to A-K stars proceeds in two steps. First, the stellar 
angular diameter <p is evaluated by comparing the IR flux observed on Earth 
in a given band to the flux predicted by a model atmosphere calculated with 
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the observed gravity and abundances and an approximate T e g (the IR flux 
does not depend sensitively on T e g). Then T e g is obtained from the total 
(integrated) flux Fboi and (j). Iteration of the procedure yields a "definite" value 
of T e g. Using IRFM, Alonso et al (1996a) derived temperatures of 475 F0-K5 
stars (T e fj in the range 4000-8000K) with internal accuracies of ~1.5%. The 
zero-point of their T e ff-scale is based on direct interferometric measures by Code 
et al (1976), and the resulting systematic uncertainty is ~1%. Accuracies of ~1% 
were obtained by Blackwell & Linas-Gray (1998), who applied IRFM to 420 A0- 
K3 stars, corrected for interstellar extinction using Hipparcos parallaxes. Both 
sets of results compare well, with differences below 0.12±1.25% for the 93 stars 
in common. 

The surface brightness method (Barnes et al 1978) was applied by Di Benedetto 
(1998) to obtain a (T e g, V-K) calibration. The calibration is based on 327 stars 
with high-precision K-magnitudes from the Infrared Space Observatory (ISO), 
Hipparcos V-magnitudes and parallaxes (the latter to correct for interstellar ex- 
tinction), and bolometric fluxes from Blackwell &: Linas-Gray (1998). First, the 
visual surface brightness Sy = V + 5 log <f> is calibrated as a function of (V-K) 
using stars with precise (ft from interferometry. Then for any star Sy is obtained 
from (V-K), (j) from Sy and V, yielding in turn T e g- from Fbd and (j). From 
the resulting (T e g, V-K) calibration, Di Benedetto derived T e g values of 537 ISO 
A-K dwarfs and giants with ±1% accuracy. The method produces results in good 
agreement with those of IRFM and is less dependent on atmosphere models. 

Multiparametric empirical calibrations of T e g as a function of the color indices 
and eventually of metallicity [Fe/H] (logarithm of the number abundances of Fe 
to H relative to the solar value) and gravity can be derived from the empirical 
determinations of the effective temperatures of the rather nearby stars. In turn, 
the effective temperature of any star lying in the (rather narrow) region of the 
H-R diagram covered by a given calibration can easily be derived (see for example 
Alonso et al 1996b). Empirical calibrations also serve to validate purely theoret- 
ical calibrations based on model atmospheres; these latter have the advantage of 
covering the entire parameter space of the H-R diagram (i.e. wide ranges of color 



indices, metallicities and gravities; see Section 3.2 later in this article). 

Spectroscopic determination of T e g is based on the analysis of chosen spectral 
lines that are sensitive to temperature; for instance the Balmer lines for stars 
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with T c gf in the interval 5000-8000 K. Because of the present high quality of 
the stellar spectra, precisions of ±50-80 K on T e ff, that correspond to the ad- 
justment of the theoretical line profile to the observed one, are commonly found 
in the literature (Cayrel de Strobel et al 1997b, Fuhrmann 1998). This sup- 
poses that theoretical profiles are very accurate, and therefore neglects the model 
atmosphere uncertainties. 

Popper (1998) used detached eclipsing binaries with rather good Hipparcos 
parallaxes, accurate radii, and measured V-flux to calibrate the radiative flux as 
a function of (B-V); he found good agreement with similar calibrations based 
on interferometric angular diameters. From the same data, Ribas et al (1998) 
derived effective temperatures (this required bolometric corrections) and found 
them to be in reasonable agreement (although systematically smaller by 2-3%) 
with T c fj derived from photometric calibrations. However the stars are rather 
distant, which implies rather significant internal errors on Mboi and T e g (a 
parallax error of 10% is alone responsible for a T e g -error of 5%). In Ribas et 
al's sample, only a few systems have a n /n <10%, and because errors on radius, 
magnitudes, and BC also intervene, only 5 systems have T e g determined to 
better than 3%. 

SURFACE GRAVITIES. If T cff and M bol are known, the radius of the star 
may be derived from the Stefan-Boltzmann law and the mass estimated from a 
grid of stellar evolutionary models, yielding in turn the value of g. This method 
has been applied to a hundred metal-poor subdwarfs and subgiants with accurate 
distances from Hipparcos (Nissen et al 1997, Fuhrmann 1998, Clementini et al 
1999). Nissen et al showed that among the various sources of errors, the error on 
distance still dominates, but pointed out that if the distance error is lower than 
20% then the error on log g may be lower than ± 0.20 dex. 

On the other hand, g can be determined from spectroscopy. Different gravities 
produce different atmospheric pressures , modifying the profiles of some spectral 
lines. Two methods have been widely used to estimate g. The first method 
is based on the analysis of the equation of ionization equilibrium of abundant 
species, iron, for instance. The iron abundance is determined from Fel lines 
that are not sensitive to gravity, and then g is adjusted so that the analysis of 
Fell lines, which are sensitive to gravity, leads to the same value of the iron 
abundance. The accuracy in log g is in the range ±0.1-0.2 dex (Axer et al 1994). 
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The second method relies on the analysis of the wings of strong lines broadened 
by collisional damping, such as Ca I (Cayrel et al 1996) or the Mg lb triplet 
(Fuhrmann et al 1997), leading to uncertainties smaller than 0.15 dex. The two 
methods often produce quite different results, with systematic differences of ~0.2- 
0.4 dex, at least when ionization equilibria are estimated from models in local 
thermodynamical equilibrium (LTE). 

Thevenin h Idiart (1999) have studied the effects of departures from LTE 
on the formation of Fel and Fell lines in stellar atmospheres, and found that 
modifications of the ionization equilibria resulted from the overionization of iron 
induced by significant UV fluxes. The nice consequence is that the gravities they 
inferred from iron ionization equilibrium for a sample of 136 stars spanning a 
large range of metallicities become very close to gravities derived either from 
pressure-broadened strong lines or through Hipparcos parallaxes. 
ABUNDANCES OF THE CHEMICAL ELEMENTS. The spectroscopic de- 
termination of abundances of chemical elements rests on the comparison of the 
outputs of model atmospheres (synthetic spectra, equivalent widths) with their 
counterpart in the observed spectra. This requires a preliminary estimate of 
T e gf and g. If high-resolution spectra are used, the line widths are very precise 
and the internal uncertainty in abundance determinations depends on uncertain- 
ties in g and T c g- , on the validity of the model atmosphere, and on the oscillator 
strengths. Error bars in the range ± 0.05-0.15 dex are typical (Cayrel de Strobel 
et al 1997b, Fuhrmann 1998). Also when different sets of [Fe/H] determinations 
are compared, the solar Fe/H ratio used as reference must be considered; values 
differing by ~0.15 dex used to be found in the literature (Axer et al 1994). This 
has resulted in long-standing difficulties in determining the solar iron abundance 
from Fel or Fell lines, because of uncertain atomic data. In a recent paper, 
Grevesse & Sauval (1999) reviewed the problem and opted for a "low" Fe- value, 
Ape = 7.50 ±0.05 (Ap e = log(np e /nH) + 12 is the logarithm of the number density 
ratio of Fe to H particles), in perfect agreement with the meteoritic value. 

Furthermore, if abundances are estimated from model atmospheres in LTE, 
perturbations of statistical equilibrium by the radiation field are neglected. Thevenin 
& Idiart (1999) found that in metal-deficient dwarfs and subgiants, the iron ove- 
rionization resulting from reinforced UV flux modifies the line widths. They 
obtained differential non-LTE/LTE abundance corrections increasing from 0.0 
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dex at [Fe/H]=0.0 to +0.3 dex at [Fe/H]= -3.0. These corrections are indeed 
supported by the agreement between spectroscopic gravities and "Hipparcos" 
gravities discussed previously. 

Helium lines do not form in the photosphere of low-temperature stars which 
precludes a direct determination of helium abundance. The calibration of the 
solar model in luminosity and radius at solar age yields the initial helium content 
of the Sun (Christensen-Dalsgaard 1982), while oscillation frequencies give access 
to the present value in the convection zone (Kosovichev et al 1992). In other 
stars, it is common to use the well-known scaling relation Y — Y p = Z^, which 
supposes that the helium abundance has grown with metallicity Z from the 
primordial value Y p to its stellar birth value Y (Y and Z represent abundances 
in mass fraction); AY/ AZ is the enrichment factor. 

a-element abundances (O, Ne, Mg, Si, S, Ar, Ca, Ti) have now been widely 
measured in metal-deficient stars. Stars with [Fe/H]<-0.5 dex generally exhibit 
an a-element enhancement with respect to the Sun ([a/Fe]) quite independent of 
their metallicity (Wheeler et al 1989, Mc William 1997). Recent determinations of 
[a/Fe] in 99 dwarfs with [Fe/H]<-0.5 from high-resolution spectra by Clementini 
et al (1999) yield [a/Fe]=+0.26±0.08 dex. 

3 RECENT THEORETICAL AND NUMERICAL PROGRESS 

Recent developments in the physical description of low and intermediate mass 
stars are briefly presented. 

3.1 Microscopic Physics 

The understanding of stellar structure benefited substantially from the complete 
re-examination of stellar opacities by two groups: the Opacity Project (OP, see 
Seaton et al 1994) and the OPAL group at Livermore (see Rogers & Iglesias 
1992). Both showed by adopting different and independent approaches, that 
improved atomic physics lead to opacities generally higher than the previously 
almost "universally" used Los Alamos opacities (Huebner et al 1977). The opacity 
enhancements reach factors of 2-3 in stellar envelopes with temperatures in the 
range 10 5 — 10 6 K. With these new opacities, (1) a number of long-standing 
problems in stellar evolution have been solved or at least lessened and (2) finer 
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tests of stellar structure could be undertaken. Since opacity is very sensitive to 
metallicity, any underlying uncertainty on metallicity may be problematic. 

Great efforts have also been invested in the derivation of low-temperature opac- 
ities, including millions of molecular and atomic lines and grain absorption that 
are fundamental for the calculation of the envelopes and atmospheres of cool stars 
(Kurucz 1991, Alexander & Ferguson 1994). 

OP and OPAL opacities have been shown to be in reasonable agreement 
(Seaton et al 1994, Iglesias & Rogers 1996); and very good agreement between 
OPAL and Alexander & Ferguson's or Kurucz's opacities is found in the domains 
where they overlap. Although some uncertainties remain that are difficult to 
quantify, the largest discrepancies between the various sets of tables do not ex- 
ceed 20% and are generally well understood (Iglesias & Rogers), making opacities 
much more reliable today than they were ten years ago. 

The re-calculation of opacities required appropriate equations of state (EOS). 
The MH&D EOS (see Mihalas et al 1988) is part of the OP project, while the 
OPAL EOS was developed at Livermore (Rogers et al 1996). In the meantime, 
another EOS was designed to interpret the first observations of very low-mass 
stars and brown dwarfs (Saumon & Chabrier 1991). OPAL and OP EOS are 
needed to satisfy the strong helioseismic constraints (Christensen-Dalsgaard & 
Dppen 1992). 

3.2 Atmospheres 

Atmospheres intervene at many levels in the analysis of observations (Section |2.2| ). 
They also provide external boundary conditions for the calculation of stellar struc- 
ture and necessary relations to transform theoretical (Mboi, T c g) H-R diagrams 
to color-magnitude (C-M) or color-color planes. Models have improved during 
the last two decades, and attention has been paid to the treatment of atomic and 
molecular line blanketing. The original programs MARCS of Gustafsson et al 
(1975) and ATLAS by Kurucz (1979) evolved toward the most recent ATLAS9 
version, appropriate for O-K stars (Kurucz 1993) and NMARCS for A-M stars 
(see Brett 1995 and Bessell et al 1998). On the other hand; very low-mass stel- 
lar atmosphere models were developed; Carbon (1979) and Allard et al (1997) 
reviewed calculation details and remaining problems (such as incomplete opacity 
data, poor treatment of convection, neglect of non-LTE effects or assumption of 
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plane-parallel geometry) . 

COLOR-MAGNITUDE TRANSFORMATIONS. Different sets of transforma- 
tions (empirical or theoretical) were used to analyze the "Hipparcos" stars. Em- 
pirical transformations have been discussed in Section |2.2| . The most recent 
theoretical transformations are compiled by Bessell et al (1998), who used syn- 
thetic spectra derived from ATLAS9 and NMARCS to produce broad-band col- 
ors and bolometric corrections for a very wide range of T e g , g and [Fe/H] values. 
These authors found fairly good agreement with empirical relations except for 
the coolest stars (M dwarfs, K-M giants). 

INTERIOR/ATMOSPHERE INTERFACE. The external boundary conditions 
for interior models are commonly obtained from T(r)-laws (r is the optical depth) 
derived either from theory or full atmosphere calculation. This method is suitable 
for low- and intermediate-mass stars (it is not valid for masses below ~0.6 M , 
Chabrier & Baraffe 1997). Morel et al (1994) and Bernkopf (1998) focused on 
the solar case where seismic constraints require a careful handling of external 
boundary conditions. Morel et al pointed out that homogeneous physics should 
be used in interior and atmosphere (opacities, EOS, treatment of convection) and 
showed that the boundary level must be set deep enough, in zones where the dif- 
fusion approximation is valid. Bernkopf discussed some difficulties in reproducing 
Balmer lines related to the convection treatment. 

3.3 Transport Processes 

CONVECTION. 3-D numerical simulations at current numerical resolution are 
able to reproduce most observational features of solar convection such as im- 
ages, spectra, and helioseismic properties (Stein & Nordlund 1998). However, 
the "connection" with a stellar evolution code is not easy, and stellar models 
still mostly rely on 1-D phenomenological descriptions such as the mixing-length 
theory of convection (MLT, Bohm-Vitense 1958). The mixing-length parameter 
omlt (ratio of the mixing-length to the pressure scale height) is calibrated so 
that the solar model yields the observed solar radius at the present solar age. 
The question of the variations of «mlt in stars of various masses, metallicities, 
and evolutionary stages remains a matter of debate (Section |]). As pointed out 
by Abbett et al (1997), the MLT can reproduce the correct entropy jump across 
the superadiabatic layer near the stellar surface, but fails to describe the detailed 



14 Yveline Lebreton 

depth structure and dynamics of convection zones. Abbett et al found that the 
solar entropy jump obtained in 3-D simulations corresponds to predictions of the 
MLT for omlt ~1.5. Ludwig et al (1999) calibrated «mlt from 2-D simulations 
of compressible convection in solar-type stars for a broad range of T e g and g- 
values. The solar omlt inferred from 3-D and 2-D simulations is close to what 
is obtained in solar model calibration. The omlt dependence with T e g and g 
of Ludwig et al can be used to constrain the range of acceptable variations of 
"mlt in stellar models (see Section 

OVERSHOOTING. Penetration of convection and mixing beyond the classical 
Schwarzschild convection cores (overshooting process) modifies the standard evo- 
lution model of stars of masses M>1.2 M Q , in particular the lifetimes (see for 
instance Maeder Ik Mermilliod 1981, Bressan et al 1981). The extent of over- 
shooting was estimated for the first time from the comparison of observed and 
theoretical MS widths of open clusters (Maeder & Mermilliod 1981), which yields 
an overshooting parameter a ov ~0.2 (ratio of overshooting distance to pressure 
scale height). As discussed in detail by Roxburgh (1997), a ov is still poorly 
constrained despite significant efforts made to establish the dependence of over- 
shooting with mass, evolutionary stage, or chemical composition (see Section Q). 
Andersen (1991) first pointed out that the simultaneous calibration of well-known 
binaries (masses and radii at 1-2%) may provide improved constraints for a ov . A 
modeling of the sample of the best-known binaries indicates a trend for a ov to in- 
crease with mass and suggests a decrease of a ov with decreasing metallicity (Ribas 
1999), although a larger sample would be desirable to confirm those trends. Fur- 
ther advances are expected from asteroseismology (Brown et al 1994, Lebreton 
et al 1995). 

DIFFUSION OF CHEMICAL ELEMENTS. Various mixing processes may occur 
in stellar radiative zones (see Pinsonneault 1997). In low-mass stars, microscopic 
diffusion due to gravitational settling carries helium and heavy elements down 
to the center and modifies the evolutionary course as well as the surface abun- 
dances. It has been proved that microscopic diffusion can explain the low helium 
abundance of the solar convective zone derived from seismology (Christensen- 
Dalsgaard et al 1993). On the other hand, turbulent mixing (resulting, for in- 
stance, from hydro dynamical instabilities related to rotation, see Zahn 1992) 
probably inhibits microscopic diffusion. Richard et al (1996) did not find any 



HIPPARCOS <fc STELLAR ASTROPHYSICS 15 

conflict between solar models including rotation-induced mixing (to account for 
Li and Be depletion at the surface) and microscopic diffusion (to account for he- 
lioseismic data). More constraints are required to clearly identify (and quantify 
the effects of) the various candidate mixing processes; this will be illustrated in 
the following sections. 

4 STUDIES OF THE BEST-KNOWN OBJECTS 

Stellar model results depend on a number of free input parameters. Some are 
observational data (mass, chemical composition and age, the latter for the Sun 
only), whereas others enter phenomenogical descriptions of poorly- known phys- 
ical processes (mixing- length parameter for convection, overshooting, etc). The 
model outputs have to be compared with the best available observational data: 
luminosity, T e g or radius, oscillation frequencies, etc. Numerous and precise ob- 
servational constraints allow assessment of the input physics or give more precise 
values of the free parameters. They may reveal the necessity to include processes 
previously neglected and in the best cases to characterize them. 

The model validation rests on (1) the nearest objects with the most accurate 
observations, (2) special objects with additional information such as stars belong- 
ing to binary systems, members of stellar clusters or stars with seismic data, and 
(3) large samples of objects giving access to statistical studies. 

4-1 Stars in Binary Systems 

Masses are available for a number of stars belonging to binary systems, allowing 
their "calibration" under the reasonable assumption that the stars have the same 
age and were born with the same chemical composition (Andersen 1991, Noels 
et al 1991). A solution is sought which reproduces the observed positions in the 
H-R diagram of both stars. Andersen (1991) claimed that the only systems able 
to really constrain the internal structure theory are those with errors lower than 
2% in mass, 1% in radius, 2% in T e g and 25% in metallicity. 

However, additional observations may sometimes cast doubts on an observed 
quantity previously determined with good internal accuracy. This occurred re- 
cently for the masses of stars in the nearest visual binary system a Centauri. 
The system has been widely modeled in the past (Noels et al 1991, Edmonds et 
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al 1992, Lydon et al 1993, Fernandes & Neuforge 1995) with the objective of get- 
ting (among others) constraints on the mixing-length parameter. At that time, 
the astrometric masses were used (internal error of 1%) but the [Fe/H]-value 
was controversial, leading to various possibilities for «mlt -values. Today the 
situation is still confused: metallicity is better assessed, but new radial velocity 
measurements yield masses higher than those derived from astrometry (by 6-7%, 
Pourbaix et al 1999). The higher masses imply a reduction in age by a factor 
of 2 and slightly different «mlt -values for the two stars. However, the orbital 
parallax corresponding to the high-mass "option" is smaller than and outside the 
error bars of both ground-based and Hipparcos parallax 7THipp- Pourbaix et al 
noted the lack of reliability of 7rHi PP given in the Hipparcos Catalogue, but since 
then it has been re-determined from intermediate data by Sderhjelm (1999) and 
is now close to (and in agreement with) the ground-based parallax. More accu- 
rate radial velocity measurements are therefore needed to assess the high-mass 
solution. 

Possible variations of ckmlt have been investigated through the simultaneous 
modeling of selected nearby visual binary systems (Fernandes et al 1998, Pourbaix 
et al 1999, Morel et al 2000). Small variations of «mlt ( n °t greater than w0.2) 
in the two components of aCen (Pourbaix et al) and i Peg (Morel et al) have 
been suggested. Fernandes et al, who calibrated 4 systems and the Sun with 
the same program and input physics, found that «mlt is almost constant for 
[Fe/H] in the range [Fe/H]0±O.3 dex and masses between 0.6 and 1.3 Mq. In 
this mass range the sensitivity of models to «mlt increases with mass (due to the 
increase with mass of the entropy jump across the superadiabatic layer) which 
makes the MS slope vary with «mlt • Also, I estimate from my models that a 
change of «mlt of ±0.15 around 1 M translates into a T e fj-change of ~40-55 K 
depending on the metallicity. On the other hand, with the solar-aMLT value 
the MS slope of field stars and Hyades stars is well fitted (Section |4.3| ). It is 
therefore reasonable to adopt the solar-aMLT value to model solar-type stars. 
For other stars, the situation is less clear. The calibration of ckmlt depends on 
the external boundary condition applied to the model, itself sensitive to the low- 
temperature opacities, and on the color transformation used for the comparison 
with observations. Chieffi et al (1995) examined the MS and red giant branch 
(RGB) in metal-deficient clusters and suggested a constancy of ckmlt from MS 
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to RGB and a decrease with decreasing Z. They found variations of «mlt with 
Z of ~ 0.2-0.4, but these are difficult to assess considering uncertainties in 
the observed and theoretical cluster sequences. On the other hand, calibration 
of ckmlt with 2-D simulations of convection gives complex results (Freytag & 
Salaris 1999; Freytag et al 1999). In particular, (1) for solar metallicity, omlt is 
found to decrease when T e g- increases above solar T e gf, and to increase slightly 
when stars move toward the RGB (by ~ 0.10-0.15) and; (2) «mlt does not vary 
importantly when metallicity decreases at solar T c g. More work is needed to go 
into finer details, and other calibrators of omlt are required, such as binary stars 
in the appropriate range of mass and with various chemical compositions. 

The modeling of a moderately large sample of binaries might give information 
on the variation of helium Y and age with metallicity Z, of great interest for 
Galactic evolution studies. The combined results for six binary systems and the 
Sun with the same program by Fernandes et al (1998) and Morel et al (2000) 
show a general trend for Y to increase with Z: Y increases from 0.25 to 0.30 
(±0.02) when Z increases from 0.007 to 0.03 (±0.002). However, the Hyades 
appear to depart from this tendency (see Section p~3| ). 

The sample of binaries with sufficiently accurate temperatures and abundances 
is still too meager to allow full characterization of physical processes. Additional 



data are needed such as observations of binaries in clusters (see Section |L3|) or 
asteroseismological measurements. 

4-2 The Nearest Disk and Halo Stars 

FINE STRUCTURE OF THE H-R DIAGRAM. Highly accurate distances for a 
rather large number of stars in the solar neighborhood were provided by Hippar- 
cos. This allowed the first studies of the fine structure of the H-R diagram and 
related metallicity effects to be undertaken. 

Among an ensemble of "Hipparcos" F-G-K stars closer than 25 pc, with error 
on parallax lower than 5%, Lebreton et al (1997b) selected stars with [Fe/H] in 
the range [—1.0, ±0.3] from detailed spectroscopic analysis (<T[f c /h] —0.10 dex, 
Cayrel de Strobel et al 1997b), Fboi and T e g from Alonso et al (1995, 1996a) with 



p b b °' ~ 2% and -tJ? 32 - ~ 1.5% (see Section 2.2 ) and not suspected to be unresolved 



binaries. Figure |l| presents the H-R diagram of the 34 selected stars: the error 
bars are the smallest obtained for stars in the solar neighborhood (aM bol are in 
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the range 0.031-0.095 with an average value (o"M bo i) —0.045 mag). The sample 
is compared with theoretical isochrones derived from standard stellar models in 
Figure |. Models cover the entire [Fe/H]-range. They account for an a-element 
enhancement [a/Fe]=+0.4 dex for [Fe/H]<-0.5 and, for non-solar [Fe/H], have a 
solar-scaled helium content (Y=Y P +Z(AY/ AZ) & ). The splitting of the sample 
into a solar metallicity sample and a moderately metal-deficient one (Figure |2|a 
and b) shows that: 

1. The slope of the MS is well reproduced with the solar cumlt , 

2. Stars of solar metallicity and close to it occupy the theoretical band cor- 
responding to their (LTE) metallicity range, while for moderately metal 
deficient stars there is a poor fit. 

In general, stars have a tendency to lie on a theoretical isochrone corresponding 
to a higher metallicity than the spectroscopic (LTE) value. This trend was already 
noticed by Axer et al (1995) but it is now even more apparent because of the 
high accuracy of the data. Helium content well below the primordial helium value 
would be required to resolve the conflict. 

This is exemplified by the star \x Cas A, the A-component of a well-known, 
moderately metal-deficient binary system that has a well-determined mass (error 
in mass of 8 per cent). The standard model (Figure ||) is more than 200 K hotter 
than the observed point and is unable to reproduce the observed T e g even if 
(reasonable) error bars are considered (Lebreton 2000). On the other hand, the 
mass-luminosity properties of the star are well reproduced if the helium abun- 
dance is chosen to be close to the primordial value, although the error bar in 
mass is somewhat too large to provide strong constraints. 

Several reasons can be invoked to explain the poor fit at low metallicities: 

1. Erroneous temperature-scale. 3-D model atmospheres could still change 
the T e gf-scale as a function of metallicity (Gustafsson 1998), but with the 
presently (1-D) available models it seems difficult to increase Alonso et al's 
(1996a) T cS by as much as 200-300 K. As noted by Nissen (1998), this scale 
is already higher than other photometric scales, by as much as 100 K. Also, 
Lebreton et al (1999) verified that spectroscopic effective temperatures lead 
to a similar misfit. 
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2. Erroneous metallicities. As discussed in Section 2.2, the [Fe/H]-values in- 
ferred from model atmosphere analysis should be corrected for non-LTE 
effects. According to Thevenin & Idiart (1999) no correction is expected at 
solar metallicity, whereas for moderately metal-deficient stars the correction 
amounts to ~0.15 dex. 

3. Inappropriate interior models. In low-mass stars, microscopic diffusion by 
gravitational settling can make helium and heavy elements sink toward the 
center, changing surface abundances as well as inner abundance profiles. 
In metal-deficient stars this process may be very efficient for three reasons: 
(1) densities at the bottom of the convection zone decrease with metallicity, 
which favors settling; (2) the thickness of the convection zones decreases with 
metallicity, making the reservoir easier to empty; and (3) metal-deficient 
generally means older, which implies more time available for diffusion. 

The two latter reasons are attractive because they qualitatively predict an in- 
creasing deviation from the standard case when metallicity decreases. As shown 
in Figure ^, a combination of microscopic diffusion effects with non-LTE Fe/H cor- 
rections could remove the discrepancy noted for //Cas A: an increase of [Fe/H] by 
0.15 dex produces a rightward shift of 80 K of the standard isochrone, represent- 
ing about one third of the discrepancy. Additionally, adding microscopic diffu- 
sion effects, according to recent calculations by Morel & Baglin (1999), provides 
a match to the observed positions. Moreover, the general agreement for solar 
metallicity stars (Figure ||a) should remain: (1) at solar metallicities non-LTE 
corrections are found to be negligible, and (2) at ages of ~5 Gyr chosen as a 
mean age for those (expectedly) younger stars, diffusion effects are estimated to 
be smaller than the error bars on T e g (Lebreton et al 1999). 

To conclude on this point, the high-level accuracy reached for a few tens of 
stars in the solar neighborhood definitely reveals imperfections in interior and 
atmosphere models. It casts doubts on abundances derived from model atmo- 
spheres in LTE, and favors models that include microscopic diffusion of helium 
and heavy elements toward the interior over standard models. Also, diffusion 
makes the surface [Fe/H]-ratio decrease by ~0.10 dex in 10 Gyr in a star like \i 
Cas (Morel & Baglin 1999), which is rather small and hidden in the observational 
error bars. In very old, very deficient stars, the [Fe/H] -decrease is expected to 
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be larger (Salaris et al 2000), which makes the relation between observed and 
initial abundances difficult to establish. In the future, progress will come from 
the study of enlarged samples reaching the same accuracies and of the acquisition 
of additional parameters to constrain the models. The knowledge of masses for 
several binaries in a narrow mass range but large metallicity range would help to 
constrain the helium abundances, while access to seismological data for at least 
one or two stars would help to better characterize mixing processes. 
STATISTICAL STUDIES. Complete H-R diagrams of stars of the solar neigh- 
borhood have been constructed by adopting different selection criteria, and have 
been compared to synthetic H-R diagrams based on theoretical evolutionary 
tracks. 

• Schroder (1998) proposed diagnostics of MS overshooting based on star 
counts in the different regions of the "Hipparcos" H-R diagram of stars 
in the solar neighborhood (d<50-100 pc). In the mass range 1.2-2 M Q , 
convective cores are small, and it is difficult to estimate the amount of over- 
shooting with isochrone shapes. Schroder suggested using the number of 
stars in the Hertzsprung gap, associated with the onset of H-shell burning, 
as an indicator of the extent of overshooting around 1.6 M ; the greater 
the overshooting on the MS, the larger the He-burning core, and in turn the 
longer the passage through the Hertzsprung gap. Actual star counts favor 
an onset of overshooting around ~1.7 M (no overshooting appears neces- 
sary below that mass), which is broadly consistent with other empirical 
calibrations (MS width, eclipsing binaries), but finer quantitative estimates 
would require more accurate observational parameters, mainly in T c g- and 
Z. 

• Jimenez et al (1998) compared the red envelope of "Hipparcos" subgiants 
(ov/V <0.15, ct(b-v) <0.02 mag) with isochrones to determine a minimum 
age of the Galactic disk of 8 Gyr, which is broadly consistent with ages ob- 
tained with other methods (white-dwarf cooling curves, radioactive dating, 
isochrones, or fits of various age-sensitive features in the H-R diagram). The 
fit is still qualitative: the metallicities of subgiants are unknown because 
of the inadequacy of model atmospheres in that region. For this reason, 
Jimenez et al investigated the isochrones in other regions, MS and clump 
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(He core burning). They calculated the variations with mass of the clump 
position for a range of metallicities in the disk, and showed that stars with 
masses from 0.8 to 1.3 MQ(ages from 2 to 16 Gyr) all occupy a well-defined 
vertical branch, the red-edge of the clump. The color of this border line is 
sensitive to metallicity, which makes it a good metallicity indicator in old 
metal-rich populations. 

• Ng & Bertelli (1998) revised the ages of stars of the solar neighborhood and 
derived corresponding age-metallicity and age-mass relations. Fuhrmann 
(1998) combined the [Mg/H]-[Fe/H] relation with age and kinematical in- 
formation to distinguish thin and thick disk stars. Several features seem to 
emerge from these studies: (1) no evident age-metallicity relation exists for 
the youngest (< 8 Gyr) thin-disk stars; some of them are rather metal-poor, 
and super metal-rich stars appear to have been formed early in the history 
of the thin disk; (2) there is an apparent lack of stars in the age-interval 
10-12 Gyr which is interpreted by Fuhrmann as a signature of the thin-disk 
formation; and (3) beyond 12 Gyr there is a slight decrease of metallicity 
with increasing age for stars of the thick disk; some of them are as old as 
halo stars. To assess these suggestions and to assist progress in the under- 
standing of the Galactic evolution scenario (see Fuhrmann 1998 for details), 
enlarged stellar samples and further improvements on age determinations 
are of course required. 

THE SUBDWARF / SUBGI ANT SEQUENCE. Hipparcos provided the very first 
high-quality parallaxes for a number of halo stars. Age determinations of the local 
halo could be undertaken, as well as comparisons with globular cluster sequences. 

Among a large sample of Population II Hipparcos halo subdwarfs, Cayrel et al 
(1997b) extracted the best-known stars with criteria similar to those adopted by 
Lebreton et al (1999) for disk stars. Stars were corrected for reddening, excluding 
stars with E(B-V)>0.05. Prior to Hipparcos, only 5 halo stars had parallax errors 
smaller than 10%; now there are 17, which represents sizeable progress. The halo 
stars are plotted in Figure |2]; subdwarfs but also subgiants are present, delineating 
an isochrone-like shape with a turn-off region. 

To make a first estimate of the age of the local halo, Cayrel et al kept 13 stars 
with the lowest error bars and spanning a narrow metallicity range ([Fe/H]=- 
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1.5±0.3), the most commonly found in the halo (Figure |5|). They found that halo 
stars, like disk stars, are colder than the theoretical isochrone corresponding to 
their metallicity. The misfit was also noted by Nissen et al (1997) and Pont et al 
(1997) in larger samples of halo stars. The discrepancy amounts to 130 to 250 K 
depending on the metallicity, and comparisons indicate that it is independent 
of the particular set of isochrones used. Again, non-LTE corrections leading to 
increased [Fe/H]-values (A[Fe/H]=+0.2 for [Fe/H]~-1.5 according to Thevenin 
& Idiart 1999), added to the effects of microscopic diffusion, can be invoked 
to reduce the misfit. Figure ||a compares Cayrel et al's sample with standard 
isochrones by PA Bergbusch & DA VandenBerg (2000, in preparation), showing 
that the subdwarf main sequence cannot be reproduced by isochrones computed 
with the LTE [Fe/H]-value, but increasing the metallicity (to mimic non-LTE 
corrections) improves the fit. Figure [|b compares the halo sequence with Profntt 
& VandenBerg's (1991) isochrones that include He sedimentation. Microscopic 
diffusion makes the isochrones redder, modifies their shape, and predicts a lower 
turn-off luminosity: the best fit with the observed sequence is achieved for an 
age smaller by 0.5-1.5 Gyr than that obtained without diffusion. Models by 
Castellani et al (1997) show that, if sedimentation of metals is also taken into 
account, including its effects on the matter opacity, the isochrone shift is smaller 
than the shift obtained with He diffusion only. 

Cayrel et al (1997b) and Pont et al (1997) estimated the local halo to be 12- 
16 Gyr old (from standard isochrones). To improve the precision more stars 
with accurate parallaxes are required. Subgiants are about 100 times rarer than 
subdwarfs, and we have only two subgiants with a n /ir <12.5% (and no subgiant 
with a n /n <5%). After Hipparcos the position of the subgiant branch is still 
poorly determined, which limits the accuracy on the age determination of the 
halo stars. 

THE ZAMS POSITIONS. The sample made of Hipparcos disk and halo stars 
spans the whole Galactic metallicity range. Figure || shows the non-evolved stars 
(M bol >5.5) of Figure g and Fi gure along with standard isochrones of various 
metallicities and solar-scaled helium ((Ay/AZ) Q =2.2). It allows a discussion of 
the position of the zero age main sequence (ZAMS) as a function of metallicity 
and implications for the unknown helium abundances. 

• MS width. Although stars generally do not lie where predicted, in particular 
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at low metallicities, the observational and theoretical MS widths are in rea- 
sonable agreement for AY/AZ =2.2. This qualitative agreement is broadly 
consistent with AY/AZ ratio of ~3±2 derived from similar measures of the 
lower MS width by Pagel & Portinari (1998) and the lower limit AY/AZ >2 
obtained by Fernandes et al (1996) from pre-Hipparcos MS. It also agrees 
with extragalactic determinations (see Izotov et al 1997) or nucleosynthetic 
predictions. 

• Helium abundance at solar metallicities. It can be noted from Figure |B] 
that there are 4 stars with Fe/H close to solar on the [Fe/H]=0.3 isochrone. 
Non-LTE [Fe/H]-corrections are negligible at solar metallicity. Microscopic 
diffusion may produce a shift in the H-R diagram: for a 0.8 Mq star of 
solar Fe/H at 5 Gyr the shift is small and comparable to the observational 
error bars (but it increases with age). These disk stars are not expected 
to be very old and the shift could instead indicate that their He-content 
is lower than the solar-scaled value. Calibration of individual objects and 
groups with metallicities close to solar indicate an increase of helium with 
metallicity corresponding to AY/AZ ~2.2 from the Sun (Lebreton et al 
1999) and AY/AZ ~2.3±1.5 from visual binaries (Fernandes et al 1998) but 
exceptions are found, such as in the (rather young) Hyades which, although 
metal-rich ([Fe/H] =0.14), appear to have a solar or even slightly sub-solar 
helium content with AY/AZ ~1.4 (Perryman et al 1998). Going into finer 
resolution would clearly require more complete data including masses for 
enlarged samples of non-evolved stars. 

• Position of metal-deficient stars. Very few metal-deficient stars have accu- 
rate positions in the non-evolved part of the H-R diagram: a gap appears for 
[Fe/H]e[-1.4,-0.3] and only 4 subdwarfs are found below [Fe/H]~-1.4. The 
empirical dependence of the ZAMS location with metallicity is impossible 
to establish for these stars, which are expected to have practically primor- 
dial helium contents. This adds to difficulties in estimating the distances of 
globular clusters (Eggen & Sandage 1962; Sandage 1970, 1983; Chaboyer et 
al 1998). 
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4-3 Stars in Open Clusters 

Hipparcos observed stars in all open clusters closer than 300 pc and in the richest 
clusters located between 300 and 500 pc providing valuable material for distance 
scaling of the Universe and for studies of kinematical and chemical evolution of 
the Galaxy. The absolute cluster sequences in the H-R diagram may be con- 
structed directly from Hipparcos distances independently of any chemical com- 
position consideration. Each sequence covers a large range of stellar masses and 
contains stars which can reasonably be considered to be born at the same time 
with similar chemical composition. Several clusters provide tests of the internal 
structure models for a wide range of initial parameters, in particular for different 
metallicities. 

THE HYADES. Obtaining high-quality astrometric data for the Hyades has been 
crucial, for it is the nearest rich star cluster, used to define absolute magnitude 
calibrations and the zero-point of the Galactic and extragalactic distance scales. 
Individual distances (mean accuracy of 5%) and proper motions were given by 
Hipparcos, providing a consistent picture of the Hyades distance, structure and 
dynamics (Perryman et al 1998). The recent determinations of the Hyades dis- 
tance modulus (m—M) are all in very good agreement while the internal accuracy 
was largely improved with Hipparcos: 

• ground-based: m — M = 3.32 ±0.06 mag (104 stars, van Altena et al 1997b) 

• HST :m-M = 3.42 ± 0.09 mag (7 stars, van Altena et al 1997a) 

• Hipparcos: m — M = 3.33 ± 0.01 mag (134 stars within 10 pc of the cluster 
center, Perryman et al 1998) 

• statistical parallaxes based on Hipparcos proper motions: m — M = 3.34 ± 
0.02 mag (43 stars, Narayanan & Gould 1999a who also showed that the 
systematic error on the parallaxes toward the Hyades is lower than 0.47 
mas). 

Greatly improved precision is seen in the H-R diagrams built with Hipparcos 
data combined with the best ground-based observations (Perryman et al 1998): 

• Figure shows 40 stars with T e g- and [Fe/H]=0.14±0.05 from detailed 
spectroscopic analysis delineating the lower part of the observational MS of 
the cluster (Cayrel de Strobel et al 1997a). 
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• Figure || is the whole H-R diagram in the (My, B-V) plane for 69 cluster 
members. Known or suspected binaries, variable stars, and rapid rotators 
have been excluded (Perryman et al 1998). Also, Dravins et al (1997) derived 
dynamical parallaxes for the Hyades members from the relation between the 
cluster space motion, the positions and the projected proper motions; these 
parallaxes are more precise (by a factor of about 2) than those directly mea- 
sured by Hipparcos, yielding in turn a remarkably well-defined MS sequence 
in the H-R diagram, narrower than that given in Figure |8] (see Figure 2 in 
Dravins et al 1997). 

• Figure |9| shows the empirical mass-luminosity (M-L) relation drawn from 
the (very accurate) masses of 5 binary systems (see caption). Nine of the 
stars are MS stars. 

Comparisons with theoretical models yield some of the cluster characteristics 
(Lebreton et al 1997a, Perryman et al 1998, Lebreton 2000): 

1. The comparison of the lowest part of the MS (Figure 0), representing the 
non-evolved stars, with theoretical ZAMS corresponding to the mean ob- 
served [Fe/H] yields the initial cluster helium content Y#=0.26±0.02 and 
metallicity Z#=0.024±0.04. Metallicity is the dominant source of the un- 
certainty on Y. 

2. The comparison of the whole observed sequence with model isochrones yields 
the cluster age. Figure ^ shows that the optimum fit is achieved with an 
isochrone of 625 ±50 Myr, Y#=0.26, Zf/=0.024 and including overshooting. 
The turn-off region (which in the Hyades corresponds to the instability 
strip of 5 Scuti stars) is rather well represented by the 625 Myr isochrone 
(see also Antonello &: Pasinetti Fracassini 1998). The quoted uncertainty 
on age only includes the contribution from visual fitting of the isochrones. 
Additional errors on age result from unrecognized binaries, rotating stars, 
color calibrations and bolometric corrections, and from theoretical models 
in particular through the parameterization of overshooting (Lebreton et al 
1995). It is therefore reasonable to give an overall age uncertainty of at least 
100 Myr. 

3. In Figure |^a the observed M-L relation is compared with the theoretical 
isochrone of 625 Myr, Yjj = 0.26, Zh = 0.024, showing an excellent agree- 
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ment. The lower part of the relation is defined by the very accurate masses 
of the two components of vB22. This system gives additional constraints 
on the Y^-value derived from ZAMS calibration. Figure |9|b illustrates how 
the positions of the two vB22 components may be used to constrain Yjj in 
the whole metallicity range allowed by observations, [Fe/H]=0.14±0.05 (see 
also Lebreton 2000). 

Furthermore, in the turn-off region of the Hyades, 5 5 Scuti stars are found 
that are quite rapid rotators (t> e sinz in the range 80-200 km.s -1 , see Antonello 
& Pasinetti Fracassini 1998). From the measurement and analysis of their oscil- 
lation frequencies and the identification of the corresponding modes by means of 
models (of same age and chemical composition), we should be able to derive the 
inner rotation profile and learn about the size of convective cores and transport 
processes at work in the interiors (Goupil et al 1996, Michel et al 1999). For 
instance, the rotation profile is related to the redistribution of angular momen- 
tum by internal motions which could be generated by meridional circulation and 
shear turbulence in a rotating medium (see Zahn 1992). On the other hand, such 
motions might induce internal mixing, and as shown by Talon et al (1997), in the 
H-R diagram rotational effects "mimic" overshooting (for instance, in a star of 
9 M©, a rotational velocity of ~ 100 km.s -1 is equivalent to an overshooting of 
«ov -0.2). 

The study and intercomparison of accurate observations of the non-pulsating 
and pulsating stars located in the instability strip should clearly provide deeper 
insight into the internal structure and properties of stars of the Hyades cluster. 
However, such analysis has to integrate the various complications related to rota- 
tion, such as the displacements in any photometric H-R diagram by amounts de- 
pending on the equatorial velocity and inclination (Maeder & Peytremann 1972, 
Perez-Hernandez et al 1999) or the splitting of oscillation frequencies, which has 
to be considered in the mode identification. 

THE PLEIADES AND OTHER OPEN CLUSTERS. The membership of stars 
in nine clusters closer than 300 pc was carefully assessed by van Leeuwen (1999a) 
and Robichon et al (1999a). Robichon et al also studied nine rich clusters within 
500 pc with more than 8 members and 32 more distant clusters. For clusters 
closer than 500 pc, the accuracy on the mean parallax is in the range 0.2-0.5 mas 
and the accuracy on the mean proper motions is of the order of 0.1 to 0.5 mas 
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per year. Results from the two groups are in very good agreement. Platais et 
al (1998) looked for (new) star clusters in Hipparcos data and found one new, a 
nearby cluster in Carina (d=132 pc) with 7 identified members. 

In order to obtain an optimal mean parallax with correct error estimates, van 
Leeuwen (1999a) and Robichon et al (1999a) worked with the Hipparcos inter- 
mediate data corresponding to each cluster, parallax and proper motion of the 
cluster center, and position of each cluster member, instead of making a straight 
average of the parallaxes of the cluster members. Stars in open clusters are lo- 
cated within a few degrees on the sky and hence were often observed in the same 
field of view of the satellite. A combined solution can be obtained from inter- 
mediate data, which allows angular correlations to be taken into account and 
the resulting parallax errors to be minimized (van Leeuwen & Evans 1998; van 
Leeuwen 1999a, b; Robichon et al 1999a). 

Mermilliod et al (1997), Robichon et al (1997), and Mermilliod (2000) com- 
pared the sequences of the various clusters in C-M diagrams derived from different 
photometric systems, and found puzzling results that are at odds with the com- 
mon idea that differences in metallicity fully explain the relative positions of the 
non-evolved parts of the MS of different clusters: 

• Some clusters have different metallicities but define the same main sequence 
in the (My, B-V) plane (Praesepe, Coma Ber, aPer, Blanco 1). For instance, 
Coma Ber has a quasi-solar metallicity while its sequence is similar to that 
of the Hyades, or of the metal-rich Praesepe. 

• Some clusters sequences (Pleiades, IC 2391 and 2602) are abnormally faint 
with respect to others, for instance Coma Ber. The metallicity of the 
Pleiades as determined from spectroscopy is almost solar, and similar to 
that of Coma Ber, but the Pleiades sequence lies (unexpectedly) ~ 0.3-0.4 
mag below the Praesepe, Coma Ber, or Hyades sequence. 

• Van Leeuwen (1999a, b) even suggested a possible (although unexpected) 
correlation between the age of a cluster and its position in the H-R diagram. 

Prior to Hipparcos, precise trigonometric parallaxes had not been obtained for 
clusters except the Hyades. Distances to open clusters were evaluated through 
the main sequence fitting technique: the non-evolved part of the (observed) clus- 
ter sequence was compared to the non-evolved part of the (absolute) lower MS 
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(ZAMS) of either (1) theoretical isochrones, (2) field stars or (3) Hyades after 
a possible correction of chemical composition differences. The magnitude differ- 
ences between absolute and apparent ZAMS directly yielded the distance modulus 
of the cluster. 

The Hipparcos distances to the 5 closest open clusters (Hyades, Pleiades, aPer, 
Praesepe, and Coma Ber) can be compared to those recently derived from MS 
fitting by Pinsonneault et al (1998); they compared theoretical isochrones, trans- 
lated into the C-M plane by means of Yale color calibrations, to observational data 
both in the (My, B-V) and (My, V-I) planes. The B-V color indice is more sen- 
sitive to metallicity than V-I (Alonso et al 1996b), so Pinsonneault et al derived 
as a by-product the value of the metallicity that gives the same distance mod- 
ulus in the two planes and compared it to spectroscopic determinations. They 
judged their distance modulii to be in good agreement with Hipparcos results 
except for the Pleiades and Coma Ber. For Coma Ber, the problem could result 
from the old VRI colors used. For the Pleiades the discrepancy with Hipparcos 
amounts to 0.24 mag, and the [Fe/H]-value derived from MS fitting in the two 
color planes agrees with the spectroscopic determination of Boesgaard & Friel 
(1990), [Fe/H]= -0.034±0.024, although values in the range -0.03 to +0.13 can 
be found in the literature. In fact, with that metallicity the Hipparcos sequence 
of the Pleiades could be reproduced by classical theoretical models, provided 
they have a high helium content. The exact value depends on the model set and 
its input physics: Pinsonneault et al found Y=0.37, Belikov et al (1998) found 
Y=0.34 but for [Fe/H]=0.10 and I find Y ~ 0.31. In any case, high helium con- 
tent is only marginally supported by observations (Nissen 1976). Pinsonneault 
et al examined other possible origins of the discrepancy (erroneous metallicity, 
age-related effects, reddening) and concluded that none of them is likely to be 
responsible for the Pleiades discrepancy. 

In parallel, Soderblom et al (1998) looked for young solar- type stars appearing 
as (anomalously) faint as the Pleiades. They found 50 field stars expected to 
be young (i.e. showing activity from Ca II H and K lines), but none of them 
lies significantly below the ZAMS. They also examined the subluminous stars 
observed by Hipparcos: they chose six stars among those lying well below the 
ZAMS, measured their spectroscopic metallicities, and found them to be metal- 
deficient with respect to the Sun with, in addition, kinematics typical of stars of 
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a thick disk or halo population. 

Soderblom et al and Pinsonneault et al concluded that, taking the Hipparcos 
results for the Pleiades at face value, it would be abnormal not to find stars sim- 
ilar to the Pleiades in the field. They inferred that the distance obtained from 
multi-color MS fitting is correct and accurate to about 0.05 mag, and concluded 
that the distance to Pleiades obtained from the analysis of Hipparcos data is 
possibly wrong at the 1 mas level, which is greater than the mean random er- 
ror. They invoked statistical correlations between right ascension and parallax 
(Pa cos s) ar i sm g from the non-uniform distribution of Hipparcos observations over 
time (and in turn along the parallactic ellipse) which affects all stars, including 
clusters. Pinsonneault et al noted that in the Pleiades the brightest stars (1) are 
highly concentrated near the cluster center and are therefore subject to spatial 
correlations which gives them nearly the same parallax, (2) have smaller a n than 
fainter stars which gives them more weight in the mean parallax, and (3) are 
those which have the highest values of cos s and also the highest parallaxes in 
the Hipparcos Catalogue. They suggest that the "true" parallax (close to that 
obtained through MS fitting) is obtained if the brightest stars with high p^ cos< 5 
are excluded from the calculation. 

Narayanan & Gould (1999b) determined the parallaxes of the Pleiades stars 
by means of Hipparcos proper motions. The resulting distance modulus has a 
rather large error bar (m — M = 5.58 ± 0.18 mag), but it is in disagreement with 
that derived directly from Hipparcos parallaxes (m — M = 5.36 ± 0.07 mag), 
and in agreement with that obtained through MS fitting (m — M = 5.60 ± 0.05 
mag). Narayanan & Gould also argue that the differences between the Hippar- 
cos trigonometric parallaxes and the parallaxes derived from Hipparcos proper 
motions reflect spatial correlations over small angular scales with an amplitude 
of up to 2 mas. 

Robichon et al (1999a, b) and van Leeuwen (1999b) have subsequently derived 
more reliable distance estimates to these clusters and performed tests that do 
not support Pinsonneault et al's conclusion. The difference between Hipparcos 
and MS fitting distance moduli is small for the Hyades (0.01 mag), whereas for 
other clusters it ranges from -0.17 mag (aPer) to +0.24 mag (Pleiades). In fact, 
except for the Hyades, the difference is always larger than the error on MS fitting 
distance modulus (0.05 mag). Robichon et al showed that while the solution 
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proposed by Pinsonneault et al improves the situation for the Pleiades, it would 
introduce new difficulties for Praesepe. By means of Monte-Carlo simulations 
of the Pleiades stars, they showed that the mean value of the Pleiades parallax 
does not depend on the correlations p^cosS- They also carefully examined distant 
stars and clusters with high p^ cosS - Through these tests, Robichon et al made the 
Hipparcos distance to Coma Ber or Pleiades more secure, and did not find any 
obvious bias on the parallax resulting from a correlation between right ascension 
and parallax, either for stars within a small angular region or for the whole sky. 

On the other hand, distances from MS fitting could be subject to higher error 
bars than quoted by Pinsonneault et al. They depend on reddening and on 
transformations from the (M^oi, T e g) to the C-M plane if theoretical ZAMS are 
used as reference (or on metallicity corrections if empirical ZAMS are compared). 
Robichon et al (1999b) compared solar ZAMS from Pinsonneault et al to those I 
calculated both in the theoretical and in the (My, B-V) planes. They showed that 
while the two ZAMS are within 0.05 mag in the theoretical H-R diagram, they 
differ by 0.15-0.20 mag in the range B-V=0.7-0.8 in the (My, B-V) plane, simply 
because different C-M transformations have been applied. Also, MS fitting often 
relies on rather old and inhomogeneous color sources (in the separate Johnson and 
Kron-Eggen RI systems) requiring transformations to put all data on the same 
(Cape-Cousins system) scale. It would therefore be worthwhile to verify the 
quality and precision of these data by making new photometric measurements of 
cluster stars. 

Let us come back to the difficult question of metallicity. As pointed out by 
Mermilliod et al (1997), photometric and spectroscopic approaches may produce 
quite different results. Metallicities have been derived recently by M Grenon 
(1998, private communication) from large sets of homogeneous observations in the 
Geneva photometric system. He obtained [Fe/H]=-0.112±0.025 for the Pleiades 
(quite different from published spectroscopic values) and [Fe/H]=0.170±0.010 
and 0.143±0.008 for Praesepe and the Hyades respectively (both in agreement 
with spectroscopy). The observed cluster sequences obtained with Hipparcos 
distances for the three clusters can be roughly reproduced by theoretical models 
computed with the photometric metallicities (and allowing for small variations 
of the helium content around the solar-scaled value) and transformed to the C-M 
plane according to the Alonso et al (1996b) and Bessell et al (1998) calibrations 
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(Robichon et al 1999b). 

In conclusion, we point out that a detailed study of the fine structure of the 
H-R diagram of the Pleiades (and other clusters) requires supplementary observa- 
tions (colors and abundances) and further progress in model atmospheres. Today 
there is no obvious solid argument against the published Hipparcos distances. In 
order to identify and understand the remaining discrepancies with stellar mod- 
els, the entire set of observed clusters has to be considered (van Leeuwen 1999a). 
Furthermore, not only the positions of the sequences in the H-R diagram but 
also the density of stars along them have to be intercompared. For instance, the 
luminosity function of young clusters exhibits a particular feature (local peak 
followed by dip) that is interpreted as a signature of pre-MS stars and might 
provide information on the initial mass function and stellar formation history 
(see the study of the Pleiades by Belikov et al 1998). On the other hand, since 
the error bars on luminosity are now small with respect to errors on color in- 
dices, stronger constraints are expected from the mass-luminosity relation, as in 
the Hyades. Observations of binaries in clusters are urgently needed and there 
is hope to detect them in the future, for, as pointed out by Soderblom et al, 
the difficult detection and measurement of visual binary orbits in the Pleiades is 
within the capabilities of experiments on board HST . 

5 RARE, FAINT, SPECIAL, OR INACCESSIBLE OBJECTS 

5.1 Globular Clusters Through Halo Stars 

Globular clusters were beyond the possibilities of Hipparcos, but the knowledge of 
distances to nearby subdwarfs gave distance estimates to a few of them through 
the MS fitting technique (Sandage 1970; Reid 1997, 1998; Gratton et al 1997; 
Pont et al 1997, 1998; Chaboyer et al 1998), comparing the non-evolved part of 
the (absolute) subdwarf main-sequence to the non-evolved part of the (observed) 
globular cluster sequence. Although simple, the technique has to be applied with 
caution: 

• Only halo subdwarfs and globular clusters with the most precise data should 
be retained. Abundances should be accurate and on a consistent scale. 
Globular cluster abundances are usually determined only for giants, while 
recent preliminary values have been obtained for subgiants in M92 (King et 
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al 1998). Abundance comparisons between (1) field and cluster stars and 
(2) dwarfs and giants have shown sometimes puzzling differences (King et 
al 1998, Reid 1999). Questions have been raised as to whether they are 
primordial or appear during evolution, but definite answers clearly require 
better spectra for all types of stars as well as spherical model atmospheres 
with better treatment of convection. The globular cluster sequence should 
be determined from good photometry well below the MS turn-off, and the 
correction for interstellar reddening has to be well estimated. Very few halo 
stars have parallaxes accurate enough to fix precisely the position of the 
ZAMS (Section 

Biases (see e.g. Lutz & Kelker 1973; Hanson 1979; Smith 1987) result- 
ing from the selection of the sample in apparent magnitude, parallax, and 
metallicity, have to be corrected for (Pont et al 1998, Gratton et al 1997); 
alternatively, samples free of biases must be selected, which implies retain- 
ing the very nearby stars with highly accurate parallaxes (Chaboyer et al 
1998, Brown et al 1997). 

Globular cluster sequence and halo sequence should (ideally) have similar 
initial chemical compositions. Because of the small number of subdwarfs 
in each interval of metallicity, it is not possible to properly establish the 
variation of the observed ZAMS position with metallicity, and to correct 
for chemical composition differences between globular clusters and subd- 
warfs empirically. Chaboyer et al (1998) found it safer to limit the method 
to globular clusters that have their equivalent in the field with the same 
[Fe/H] and [a/Fe] content. Gratton et al (1997) and Pont et al (1998) 
applied theoretical color corrections to the subdwarf data to account for 
metallicity differences with globulars. In addition, element sedimentation 



might introduce further difficulties, as already mentioned in Section 4J2. As 
pointed out by Salaris & Weiss (2000), the present surface chemical com- 
position of field subdwarfs no more reflects the initial one if microscopic 
diffusion has been efficient during evolution, while, in globular cluster gi- 
ants which have undergone the first dredge-up, the chemical abundances 
have been almost restored to the initial ones. 

Unresolved known or suspected binaries can introduce errors in the defi- 
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nition of the ZAMS position. Chaboyer et al and Gratton et al excluded 
them whereas Pont et al applied an average correction of 0.375 mag on their 
position, a procedure that has been criticized by Chaboyer et al and Reid 
(1998). 

• Evolved stars have to be excluded (there is no certainty that globular clusters 
and halo dwarfs have exactly the same age). From theoretical models it is 
estimated that stars fainter than My— 5.5 are essentially unevolved. 

The number of globular clusters studied by the different authors varies because 
of the different criteria and techniques chosen to select the subdwarfs samples. 
Nevertheless they all agree on the general conclusion that globular cluster dis- 
tances derived from MS fitting are larger by ~5-7% than was previously found. 
Chaboyer (1998) calculated an average of distances to globular clusters obtained 
with different methods (MS fitting, astrometry, white dwarf sequence fitting, cali- 
bration of the mean magnitudes of RR Lyrae stars in the Large Magellanic Cloud, 
comparison with theoretical models of horizontal branch stars, statistical parallax 
absolute magnitude determinations of field RR Lyrae from the Hipparcos proper 
motions, etc) and noted that the distance scale is larger (by 0.1 mag) than his 
pre-Hipparcos reference. 

It is worth pointing out that the statistical parallax method alone favors a 
shorter distance scale (by ss0.3 mag with respect to MS fitting result). As re- 
viewed by Layden (1998) and by Reid (1999), the statistical parallax method 
was applied by independent groups who found concordant results. However, the 
absolute magnitudes My(RR) of the halo RR Lyrae derived from the statistical 
parallax method (on the basis of Hipparcos proper motions and of radial veloci- 
ties) are also ss 0.3 mag fainter than the magnitudes obtained through a method 
directly based on Hipparcos parallaxes (Groenewegen & Salaris 1999); these lat- 
ter being in turn in good agreement with the MS fitting result. As discussed 
thoroughly by Reid (1999), there are several difficulties related to the My(RR)- 
calibration and to its comparison with other distance calibrations that still hinder 
the coherent and homogeneous understanding of the local distance scale. 

In Caputo's (1998) and Reid's (1999) reviews the ages of globulars are dis- 
cussed. After Hipparcos, ages of globular clusters are reduced by typically 2-3 
Gyr, because of both larger distances and improvements in the physics of the 
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models, mainly in the equation of state and in the consideration of microscopic 
diffusion. Present ages are now in the range 10-13 Gyr, which can be compared 
to the previous interval of 13-18 Gyr (see VandenBerg et al's review, 1996). 

Chaboyer et al (1998) claimed that the (theoretical) absolute magnitude My and 
lifetimes of stars at the MS turn-off (T-O) in globular clusters are now well un- 
derstood, since the physics involved is very similar to that of the Sun, which 
is in turn well constrained by seismology. In particular, My (T-O) is quite in- 
sensitive to uncertainties related to model atmospheres or convection modeling 
(see also Freytag Sz Salaris 1999). Chaboyer et al suggested that no significant 
changes (more than ~5%) in the derived ages of globular clusters are expected 
from future improvements in stellar models. Conversely, distance and abundance 
determinations are far from definite, and the quasi-verticalness of isochrones in 
the T-0 region makes the determination of My (T-O) difficult (see Vandenberg 
et al 1996). Further revision of the ages is therefore not excluded. Also a global 
agreement between an entire globular cluster sequence and the corresponding 
model isochrone is far from being reached. 

Age of the Universe. The ages of the oldest objects in the Galaxy, the most 
metal-poor halo or globular cluster stars, provide a minimum value for the age 
of the Universe Ttj. Globular cluster ages (10-13 Gyr) from comparisons of 
isochrones with observed My(T-O) are presently the most reliable, but two in- 
dependent methods look promising: 

• Thorium (half- life 14.05 Gyr) has been detected and measured by Sneden 
et al (1996) in an ultra-metal-poor giant, too faint for observation by Hip- 
parcos, and the star radioactive decay age is estimated to be 15±4 Gyr. 
In the future, such observations of more stars and the possible detection of 
Rhenium and Uranium could provide strong constraints for Ttj. 

• Observations of (faint) white dwarfs (WD) in globular clusters are now 
within reach of experiments on board HST , and a lower limit to the age of 
WD in M4 of ~9 Gyr has been derived from a comparison with theoretical 
WD cooling curves (Richer et al 1997). Future access to cooler and fainter 
objects will better constrain Ttj. 

According to Sandage & Tammann (1997) and Saha et al (1999) the Hubble 
constant Hq should be in the range 55 ± 5km.s _1 .Mpc _1 , which implies Ty = 
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^Hq 1 ~ 11 — 13.5 Gyr, indicating that no strong discrepancy with the age of the 
oldest known stars remains. 

5.2 Variable Stars 

I shall not discuss the revisions of the distance scale based on pulsating stars (RR 
Lyrae, Cepheids, Miras, high-amplitude 5 Scuti stars) because this topic has been 
extensively reviewed by Caputo (1998) and Reid (1999). 

Both new insight as well as new questions about the physics governing pulsat- 
ing stars have been generated from the combination of Hipparcos distances with 
asteroseismic data. When the magnitude of a star is modified and its error box 
reduced, the mass and evolutionary stage attributed to the star may be modified. 
For variable stars, a different evolutionary stage may give a drastically different 
eigenmode spectrum, and in turn may change the mode identification and aster- 
oseismic analysis (see Liu et al 1997). H0g & Petersen (1997) showed that for 
the two double-mode, high-amplitude 5 Scuti variables SX Phe and AI Vel, the 
masses derived on one hand from stellar envelope models and pulsation theory, 
and on the other hand from the position in the H-R diagram through stellar 
evolution models are in nice agreement if the Hipparcos parallaxes (accurate to 
5-6%) are used. Further implications of Hipparcos distances on the understand- 
ing of S Scuti stars, A Bootis, and rapidly oscillating Ap stars have been discussed 
in several papers (see for instance Audard et al 1998, Viskum et al 1998, Paun- 
zen et al 1998, Matthews et al 1999), whereas the physical processes relevant to 
the Asymptotic Giant Branch and pulsation modeling of Miras and Long Period 
Variables were examined by Barthes (1998, see also references therein). 

5.3 White Dwarfs 

The white dwarf (WD) mass-radius (M-R) relation was first derived by Chan- 
drasekhar (1931) from the theory of stars supported by the fully degenerate 
electron gas pressure. It has been refined by Hamada & Salpeter (1961), who 
calculated zero-temperature (fully degenerate) WD models of different chemical 
composition (He, C, Mg, Si, S, and Fe) and by Wood (1995), who calculated WD 
models with carbon cores and different configurations of hydrogen and/or helium 
layers and followed the thermal evolution of WD as they cool. Although theoreti- 
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cal support is strong, it has long been difficult to confirm the relation empirically 
because of (1) the very few available WD with measures of masses and radii, (2) 
the size of the error bars and (3) the intrinsic mass distribution of the WD, which 
concentrates them in only a small interval around 0.6 M (Schmidt 1996). 

The M-R relation, assuming that WD have a carbon core, is a basic underlying 
assumption in most studies of WD properties. It serves to determine the mass of 
WD, and in turn their mass distribution and luminosity function. It is important 
because WD feature in many astrophysical applications such as the calibration 
of distances to globular clusters (Renzini et al 1996) or the estimate of the age of 
Galactic disk and halo by means of WD cooling sequences (see Winget et al 1987, 
d'Antona Sz Mazzitelli 1990). The more precisely the M-R relation is defined by 
observations, the better the tests of theoretical models of WD interiors that can 
be undertaken. These include tests of the inner chemical composition of WD, 
thickness of the hydrogen envelope of DA WD, or the characterization of their 
strong inner magnetic fields. 

Depending on the white dwarf considered, the empirical M-R relation can be 
obtained by different means: 

1. Surface brightness method. If T c g and log g are determined (generally 
from spectroscopy), then model atmospheres allow calculation of the energy 
flux at the surface of the star, which when compared to the flux on Earth 
yields the angular diameter cf> (see Schmidt 1996). The radius R is obtained 
from the parallax and cf>, M is deduced from R and log g. This method 
requires high-resolution spectra and largely depends on model atmospheres. 

2. Gravitational redshift. The strong gravitational field at the surface of a 
WD causes a redshift of the spectral lines, the size of which depends on the 
gravitational velocity t> grs = (c is the speed of light). If be 
measured and the gravity is known, then M and R can easily be obtained 
independently of the parallax. v grs can only be measured in WD members 
of binary systems, common proper motion pairs (CPM), or clusters because 
the radial velocity is required to distinguish the gravitational redshift from 
the line shift due to Doppler effect. Also, very high-resolution spectra are 
needed. 

3. WD in visual binary systems. Masses may be derived directly from the 
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orbital parameters through the Kepler's third law, provided the parallax is 
known. Radii are derived from the knowledge of T e g and distances. 

More than 15 years ago, when the Hipparcos project began, uncertainties on 
WD ground-based parallaxes were at least 10 mas. During the last 10 years, 
due to great instrumental progress, parallax determinations were improved by a 
factor of about 2, and more accurate atmospheric parameters T e g, g and t> grs were 
obtained. In the meantime, Hipparcos observed 22 white dwarfs (11 field WD, 4 
WD in visual binaries, and 7 in CPM systems) among which 17 are of spectral 
type DA. Although they are close to the faint magnitude limit of Hipparcos, the 
mean accuracy on their parallaxes is a n ~3.6 mas (Vauclair et al 1997). 

Vauclair et al (1997) and Provencal et al (1998) studied the whole sample of 
WD observed by Hipparcos. The M-R relation is narrower and most points are 
within 1 a of Wood's (1995) evolutionary models of WD with carbon cores and 
hydrogen surface layers. The theoretical shape is still difficult to confirm because 
of the lack of objects in the regions of either high or low mass. Furthermore, 
the error bars are still too large to distinguish fine features of the theoretical 
models, such as between evolutionary and zero-temperature sequences or thick 
and thin hydrogen envelopes (Vauclair et al 1997), except for some particular stars 
(Shipman et al 1997, Provencal et al 1998). Other effects such as alterations due 
to strong internal magnetic fields are not yet testable (Suh & Mathews 2000). 

WD in binary systems. Prior to Hipparcos, Sirius B was the only star roughly 
located on the expected theoretical M-R relations; the others (Procyon B, 40 
Eridani B and Stein 2051) were at least 1.5cr below the theoretical position (see 
Figure 1 in Provencal et al 1998). After Hipparcos, as shown by Provencal et al, 
the error on the radius is dominated by errors on flux and T e g. On the other 
hand, the parallax error still dominates the error on mass, except for Procyon 
where the error on the component separation plays a major role. 

• Sirius B is more precisely located on a Wood's (1995) M-R relation for a 
DA white dwarf of the observed T c g with a thick H layer and carbon core 
(Holberg et al 1998). Also compatible with Wood's thick H layer models 
is the position of V471 Tau, a member of an eclipsing binary system for 
which the Hipparcos parallax supports the view that it is a member of the 
Hyades (Werner & Rauch 1997, Barstow et al 1997). The mass of 40 Eri B 
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increased by 14%. The star is now back on Hamada & Salpeter's (1961) M- 
R relation for carbon cores, making it compatible with single star evolution 
(Figure 1 by Shipman et al 1997) and it does not appear to have a thick 
H layer. Sirius B (the most massive known WD) and 40 Eri B (one of the 
less massive) nicely anchor the high-mass and low-mass limits of the M-R 
relation. 

• The case of Procyon B remains puzzling. The position is not compatible 
with models with carbon cores, and would be better accounted for with 
iron or iron-rich core models (Provencal et al 1997). Provencal et al (1998) 
also examined seven white dwarf members of CPM pairs (more distant and 
fainter than WD in visual binaries) with Hipparcos distances and gravita- 
tional redshift measurements. They showed that two of them also lie on 
theoretical M-R relations corresponding to iron cores. This is not predicted 
by current stellar evolution theory, and further work is required to clarify 
this problem. 

In conclusion, better distances from Hipparcos and high-resolution spectroscopy 
has allowed better assessment of the theoretical M-R relation for white dwarfs, 
and has shown evidence for difficulties for a few objects that do not appear to 
have carbon cores. Future progress will come from further parallax improvements 
and from better T e g , gnitudes, and orbital parameters for visual binaries. 

A better understanding of the atmospheres is required: convection plays a role 
in the cooler WD, additional pressure effects due to undetectable helium affect 
the gravity determination, and incorrect H layer thickness estimates change the 
mass attributed to WD. Further information coming from asteroseismology or 
spectroscopy would help. 

6 FUTURE INVESTIGATIONS 

Hipparcos has greatly enlarged the available stellar samples with accurate and 
homogeneous astrometric and photometric data. To fully exploit this new in- 
formation, many studies have been undertaken (several hundred papers devoted 
to stellar studies based on or mentioning Hipparcos data can be found in the 
literature), therefore the present review could not be fully exhaustive. 

The Hipparcos mission succeeded in clarifying our knowledge of nearby objects, 
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and allowed first promising studies of rarer or farther objects. After Hipparcos, 
the theory of stellar structure and evolution is further anchored, and some of its 
physical aspects have been better characterized. For instance, new indications 
that the evolution of low-mass stars is significantly modified by microscopic dif- 
fusion have been provided by fine studies of the H-R diagram, and consequences 
for age estimates or surface abundance alterations have been further investigated. 
On the other hand, Hipparcos left us with intriguing results that raised new ques- 
tions. For example, the unexpected position of the white dwarf Procyon B on 
a theoretical mass-radius relation corresponding to iron cores is still not under- 
stood. 

Today, uncertainties on distances of nearby stars have been reduced signif- 
icantly such that other error sources emerge to dominate, hindering further 
progress in the fine characterization of stellar structure. Progress on atmosphere 
modeling is worth being pursued, for it has implications for observational param- 
eters (effective temperatures, gravities, abundances, bolometric corrections), the- 
oretical models (outer boundary conditions), and color calibrations. A thorough 
theoretical description of transport processes (convection, diffusion) and related 
effects (rotation, magnetic fields) is needed to improve stellar models, as well 
as further improvements or refinements in microscopic physics (low-temperature 
opacities, nuclear reaction rates in advanced evolutionary stages). 

What is now needed from the observational side is (1) enlarged samples of rare 
objects (distant objects, faint objects or objects undergoing rapid evolutionary 
phases), (2) an increased number of more "common" objects with extremely 
accurate data (including masses), and (3) a census over all stellar populations. 

These goals should be (at least partially) achieved by future astrometric mis- 
sions. The NASA Space Interferometry Mission (SIM), scheduled for launch in 
2006, will have the capability to measure parallaxes to 4 ^arcsecond and proper 
motions to 1-2 ^arcsecond per year down to the 20 th magnitude which repre- 
sents a gain of three orders of magnitude with respect to Hipparcos (Peterson & 
Shao 1997). The ESA candidate mission GAIA is dedicated to the observation 
of about one billion objects down to V~20 mag (and typical a w ~10 /uarcsecond 
at V=15 mag). GAIA will also provide multi-color, multi-epoch photometry for 
each object, and will give access to stars of various distant regions of the Galaxy 
(halo, bulge, thin and thick disk, spiral arms). It is aimed to be launched in 2009 
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(Perryman et al 1997b). 

Asteroseismology has already proved to be a unique tool to probe stellar in- 
teriors. Space experiments are under study. The first step will be the COROT 
mission (aimed to be launched in 2003), designed to detect and characterize oscil- 
lation modes in a few hundred stars, including solar-type stars and 5 Scuti stars 
(see Baglin 1998). 
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Figure 1: Hipparcos H-R diagram of the 34 best-known nearby stars. The paral- 
lax accuracies a^/n are in the range 0.003-0.041. Bolometric fluxes and effective 
temperatures are available from Alonso et al's (1995, 1996a) works ( p bo1 ~ 2% 
and ~ 1.5%, see Section ^L% ). Resulting (XM bol are in the range 0.031-0.095 



mag (from Lebreton et al 1999). 
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Figure 2: The sample of Figure |T| split into two metallicity domains. Figure ^a 
shows stars with [Fe/H] close to solar ([Fe/H] e[-0.45, +0.25]). Figure shows 
moderately metal-deficient stars ([Fe/H] £ [-1.05,-0.45]). Theoretical isochrones 
are overlaid on the observational data. Figure |^a: the lower isochrone (10 Gyr) 
is for [Fe/H]=-0.5, Y=0.256 and [a/Fe]=0.4; the upper isochrone (8 Gyr) is 
for [Fe/H]=+0.3, 1==0.32 and [a/Fe]=0.0; the dashed line is a solar ZAMS 
(«mlt =1-65, Y"0=O.266 and Z =O.O175). The brightest star is the young star 7 
Lep. Figure ||b: 2 isochrones (10 Gyr) with [a/Fe]=0.4, the lower is for [Fe/H]=- 
1.0, y=0.236 and the upper for [Fe/H] =-0.5, y=0.256. All stars but one are 
sitting above the region defined by the isochrones. (AY/AZ) & =2.2 is obtained 
with Balbes et al's (1993) primordial helium 1^=0.227 (from Lebreton et al 1999). 
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Figure 3: H-R diagram for the unevolved moderately metal deficient stars of 
Figure ^ (mean LTE metallicity [Fe/H]LTE = —0.72, mean non-LTE value 
[Fe /H] nlte = —0.57, see text). Full and dashed lines are standard isochrones (10 
Gyr) computed with, respectively, the [Fe/H]LTE and [Fe/H]NLTE values. The 
dot-dashed isochrone (10 Gyr) includes He and heavy elements sedimentation: 
at the surface it has [Fe/H]NLTE = —0.57 but the initial [Fe/H] was ss-0.5. (from 
Lebreton et al 1999). 
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Figure 4: Hipparcos H-R diagram of the 32 halo stars with a n /ir <0.22 (the 
parallax accuracies a^/n are in the range 0.007-0.214). Bolometric fluxes and 
effective temperatures are available from Alonso et al's (1995, 1996a) works (see 
the caption for Figure ||). Resulting crM bol are in the range 0.03-0.48 mag. 
A bunch of subgiants emerges with an isochrone-like shape (from Cayrel et al 
1997b). 
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Figure 5: Hipparcos H-R diagram of 13 halo stars with [Fe/H]LTE = —1-5 ± 0.3 
and ct^/tt <0.12 (the parallax accuracies cr n /ir are in the range 0.01-0.12). Bolo- 
metric fluxes and effective temperatures are available from Alonso et al's (1995, 
1996a) works (see the caption for Figure |l|). Resulting <JM bol are in the range 
0.03-0.26 mag. All isochrones were kindly provided by DA Vandenberg. Fig |5|a 
illustrates the effect of a non-LTE correction of +0.2 dex in [Fe/H] as inferred 
from PA Bergbusch & DA VandenBerg's (2000, in preparation) models (see text): 
the dashed line is a standard isochrone of 12 Gyr ([a/Fe]=+0.3, y~0.24) with 
[Fe/H] =-1.54 (LTE value), and the full lines are isochrones with [Fe/H]=-1.31 
(non-LTE value) of 12 Gyr (upper line) and 14 Gyr (lower line). Fig ||b il- 
lustrates the effect of microscopic diffusion of He as inferred from Proffitt's & 
VandenBerg's (1991) models: isochrones ([Fe/H]=-1.3 and [O/Fe]=0.55), of 12 
Gyr with (full line) and without (dashed-line; upper line 12 Gyr, lower 14 Gyr) 
microscopic diffusion are plotted. 
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Figure 6: "Hipparcos" H-R diagram of non-evolved stars with a^/ir <0.10. Each 
star is labeled with its [Fe/H]-value. Standard isochrones are plotted with, from 
left to right, [Fe/H]=-2.0, -1.5, -1.0, -0.5, 0.0, 0.3 (from Lebreton 2000). 
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Figure 7: Hipparcos H-R diagram for 40 selected low MS stars in the Hyades. 
The mean [Fe/H] is 0.14± 0.05. Stars with error bars are not suspected to be 
double or variable. Internal errors on T e g are in the range 50-75 K. Double 
or variable stars are also indicated: objects resolved by Hipparcos or known 
to be double systems are shown as circles, triangles denote objects with either 
detected photocentric acceleration or objects possibly resolved in photometry, and 
V means spectroscopic binary or radial velocity variable. Theoretical ZAMS loci 
are given for the Hyades (dashed line, Y=0.26 Z=0.024) and solar (dotted line, 
y=0.266 Z=0.0175) chemical compositions (from Perryman et al 1998). 
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Figure 8: Hipparcos C-M diagram of the Hyades. V and B-V values are from 
the Hipparcos catalogue (o-(b-v) < 0.05 mag). The loci of ZAMS and theoretical 
isochrones calculated with overshooting (a ov =0.2) are indicated. Arrows indicate 
the position of the components of the binary system 9 2 Tau used for the age 
determination (from Perryman et al 1998). 
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Figure 9: Figure |9|a: the Hyades empirical and theoretical mass- luminosity rela- 
tions. Masses are from Peterson &; Solensky (1988, vB22), Torres et al (1997a, b, 
c, Finsen 342, 8 1 , 6 2 Tau) and Sderhjelm (1999, 51 Tauri from Hipparcos data). 
The isochrone is for 625 Myr, taken from the same calculations used in Figure ||. 
Figure |9|b illustrates how the precise positions of the members of vB22 allow to 
discriminate between different (Y, [Fe/H]) values. 



